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Abstract 

We formulate a systematic Bethe-Ansatz approach for computing bound-state ("breather"' 
S matrices for integrable quantum spin chains. We use this approach to calculate the 
breather boundary S matrix for the open XXZ spin chain with diagonal boundary 
fields. We also compute the soliton boundary S matrix in the critical regime. 



1 Introduction 



A common feature of integrable models is the existence of bound states for a certain range of 
the coupling constant. A well-known example is the sine- Gordon/massive Thirring model, 
which in the attractive regime (0 < f3 2 < Att) exhibits soliton-antisoliton bound states called 
"breathers." (See e.g. [jlj and references therein.) The direct Bethe-Ansatz calculation of 
exact scattering matrices for both solitons (also known as "kinks" or "holes" ) and breathers 
was pioneered by Korepin 0. Andrei and Destri || later systematized such S matrix 
calculations for the solitons. We develop here a corresponding systematic approach for 
computing S matrices for the breathers. In particular, we give a direct calculation of the 
breather boundary S matrix for the open XXZ spin chain with diagonal boundary fields 
Our results coincide with the bootstrap results for the boundary sine-Gordon model 
H with "fixed" boundary conditions which were obtained by Ghoshal [fj]]. We also give 
a direct computation of the soliton boundary S matrix in the critical regime U,|8j using 
the method developed in P|,[]10|. Although we focus on the XXZ chain, we expect that 
our method of computing bound-state S matrices should be applicable to other integrable 
quantum spin chains. 

Bulk calculations are generally more straightforward than corresponding boundary cal- 
culations. We therefore first formulate in Section 2 the method of computing breather S 
matrices for the case of bulk (two-particle) scattering in the closed XXZ chain, and thereby 
reproduce the well-known results [|I|], @, [fElH . In Section 3 we turn to the open XXZ chain. 
We compute the breather boundary S matrix, and find agreement with the bootstrap results 
provided a certain identification of boundary parameters is made. In order to further check 
this identification, we also compute the soliton boundary S matrix. A brief comparison of 
our approach with that of other authors is given in Section 4. 



2 Closed XXZ chain 



In this Section we consider the periodic anisotropic Heisenberg (or "closed XXZ" ) spin chain 
in the critical regime, whose Hamiltonian is given by ||ll||,||12|,fil 



e N 

n = jJ2 {«+x + «+i + cos ^ - l ) } > °N+i = °i , (2.1; 



n=l 



with < fi < | and e = ±1. We also assume that the number of spins (N) is even. It can be 



shown (see e.g. [|5|,||ll||,[nj|) that the kink S matrix coincides with the sine-Gordon soliton 
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S matrix [jl], provided that the sine-Gordon coupling constant [3 2 is identified as 

„ 9 } 8 (it — u) for e = +1 . 

P s f 1 • (2 - 2) 

I o// tor e = — 1 

Since we restrict the anisotropy parameter /i to the range (0,|), it follows that the case 
e = +1 corresponds to the "repulsive" regime (in < f3 2 < 8n) of the sine-Gordon model 
in which there are no bound states, while e = — 1 corresponds to the "attractive" regime 
(0 < f3 2 < 47r) in which bound states do exist. 

Choosing the pseudovacuum to be the ferromagnetic state with all spins up, the algebraic 
Bethe Ansatz [[HJ can be used to construct simultaneous eigenstates of the Hamiltonian, 
momentum, and S z . The corresponding eigenvalues are given by [] 

M j 

E = -esinVy — , /0 x , , (2.3) 

^[ cosh{2jj,A a ) — cos fi 

e m S mh fi f\ + i) 

P = 7rM0(-e) + T ^ log 7 ~T (mod 2tt) , (2.4) 

i a =i sinh/i (A a - ~J 

S z = f-M, (2.5) 
where {Ai , . . . , Am} are solutions of the Bethe Ansatz equations 

M 

ei(\ a ;iJ,) N = Y[ e 2 (A a - A/? ;//) , a = l,---,M, (2.6) 



/3=1 
I3^a 



where 



sinh/ifA + fl 

e n (A;//) = ) (2.7) 

sinh/i(A-f) 

Moreover, 6? (a;) is the Heaviside unit step function. 

For the analysis that follows, it is convenient to also introduce the following notations: 

m cosh/xfA + f) 

g n (X ; fi) = e n (A ± — ; fi) = - \ 2 J , (2.8) 
2/i cosh fi f A - y J 

g„(A;/i) = 7r-Mloge„(A;/i) , 

r n (A;/i) = ilogg n (X;fi) , (2.9) 



lr rhe dependence on e in the formulas that follow is explained in Appendix A. 
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1 d fi sin(n/z) 

a n (A;/i) = — — q n (A;n)- 



2tt dX ' 7r cosh(2/iA) — cos(n/i) ' 
7 / \ \ Id / , > n sin(n/i) 

& nA;/i = — — r n {X;fj,) = — — -. 2.10) 

Z7T d\ 7T cosn(2/iA) + cos(nj2) 

The latter functions, which have the periodicity n — > n + — , have the following Fourier 
transforms Q: 

sinh ( (i/ — n)% ) 

a n (w;^) = v , r y , 0<n<2z/, (2.11) 

sinh (f ) 

sinh f ^ 

b n (u;n) = — : )-=-(, 0<n<z/, (2.12) 



sinh (if) 

sinh((n-2i/)| 
sinh (J ) 



u<n<2u, (2.13) 



where u — - > 2. 



2.1 Ground state 

In order to study the breathers, we must consider the attractive case e = — 1. The ground 
state lies in the sector with N even, and is characterized by a "sea" of M = y negative-parity 
1-strings (i.e., roots of the form A + y 1 , where the "center" A is real) ]13[. We briefly review 
the procedure for determining the root density, which describes the distribution of roots in 
the thermodynamic (N — > oo) limit. The Bethe Ansatz Eqs. ( |2.6| ) for the ground state are 

M 

gi (\ a ;fx) N = J] e 2 (A a -A^;/i), a = l,---,M, (2.14) 

(3=1 

with {A Q } real. By taking logarithms, these equations can be rewritten as 

h(\ a ) = J a , a = 1, -",M, (2.15) 

2 Our conventions are 

/oo -i rOO 

e- A /(A) dA , /(A) = — / A />) , 

and we use * to denote the convolution 

/oo 
/(A - A') g(X') d\' . 
-OO 



3 



where the so-called counting function h(X) is given by 



MA) 



1 

2^ 



M 



Nn(X ; fj) - 92(A - Ays ;/*) 

/3=1 



(2.16) 



and { J a} are certain integers or half- integers. The sign of the counting function is chosen so 
as to make it a monotonically increasing function of A. The root density er(A) is defined by 

1 d 



(2.17) 



so that the number of X a in the interval [A , A + dX] is Na(X)dX. It is a positive function 
by virtue of the monotonicity of the counting function. Passing from the sum in h(X) to an 
integral, we obtain a linear integral equation for the root density 



a(X) = -h(X ]fi)+ dX'a(X') a 2 (X - X' ; //) . 



(2.18) 

Solving this equation by Fourier transforms using Eqs. ( |2.11|) , (|2.12|) , we conclude that the 
root density for the ground state is given by 

1 



* A = * A = 7T- / 

Z7T J-oo 



dtu e 



-iu>X 



situ 



2(1/ - 1) cosh 



(2.19) 



where 



-l + a 3 (u;n) 2cosh((i/-l)|) ' 



(2.20) 



We verify the consistency of this procedure by computing the value of M from the root 
density: 

M=^1 = N dA a(A) = iVs(O) = — , (2.21) 

.. i J —oo 2 



a=l 



and hence, the state indeed has S z = 0. The energy and momentum are 

„ 7TSinU JH, 27T 7TSinU JH_ . 7T sin Lit roc 

£ 9r = ^]T ai (A Q + — ;/i) = -£&i(A a ;^) = ^JV/ cZA s(A) b x {X ; //) 



sin // 
4// 



AT 



da; 



smh (|) 



cosh((i/- l)f) sinh (*f*) 



a=l 



/7T 



<?i(A a + — ;//) - vr 



a=l 



TriV 



(mod 27r) . 



/oo 
dA s(A) ri(A;^) 
-oo 

(2.22) 
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2.2 Two-breather state 



As for the massive Thirring/ sine-Gordon model , the XXZ chain in the attractive regime 
has two classes of excitations above the ground-state sea: holes which correspond to solitons, 
and strings which correspond to soliton-antisoliton bound states, i.e., breathers. The n th 
breather corresponds to a positive-parity ra-string; i.e., a set of n roots of the Bethe Ansatz 
Eqs. of the form 

\M = \ Q +U n + l-2l) , l = l,...,n, (2.23) 

where the center A is real. In particular, the fundamental breather (n = 1) corresponds to a 
real root of the Bethe Ansatz equations. Breather states exist only for n G {1 , . . . , [v] — 1}, 
where [x] denotes integer part of x (see P|,fT3|j). 

We consider now an excited state consisting of two breathers Xi , X^ 2 ' 1 ^ (with centers 
X\ and AS], respectively) in the sea, again with N even. The Bethe Ansatz Eqs. ( B.6|) now 
imply 

M~ 2 rip 

9i(\ a ;») N = -IW A "-W)n n^(A Q -A^ ) ;/i), a = l,---,Mr, 

(3=1 f3=ll =l 

(2.24) 



Mi ri2 



f3=l h=l 

(2.25) 

where Mf is the number of roots in the sea, and Ai , . . . , X M - are real. 

The first set of equations ( 2.24 ), which describes the (distorted) sea, implies the counting 
function 

1 I Ml_ 2 n f3 , ( ^ 

} • (2-26) 



MA) = - r pr!(A ; ai) - £fc(A - A, ; 0) - £ £ r 2 (A - A< n <" W ; /x) 



The corresponding root density (|2.17| ) is therefore given by 



i 2 "3 

<r(A) = s(A) - ± £ £ i^(A - X^ M ) (2.27) 
= aW-^E^CA-Aj), (2.28) 

iV /3=1 
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where the Fourier transform of K n (X) is given by 

K-ifa ; fi) + k + i(u; ■ /x) sinh (?f) cosh (f ) 



-l + a 2 (u;;/x) sinh (f ) cosh 



(2.29) 



keeping in mind that n < v — 1. A calculation analogous to ( 2.21 ) shows that = 
y — rii — rx 2 , and therefore, the breathers have S" 2 = 0. The energy of the state is given by 



E 



7T sin /X 

7rsin /x 
H 

E gr + 



E«i(A a + ^iM) + EE a x {\^ ; /i) 



2 n D 



Q = l 



iV / dAa(A) &i(A;//) + $> na (A°;//) 

. a=l 



7T sin /X 



EZ > 

where the Fourier transform of e n (\) is given by 
e n (u) = a n (u ; /x) - lf n (a>)&i(a; ; /x) 



(2.30) 



cosh ( (v — n — 1) 



cosh ((u — 1] 



(2.31) 



which is invariant under n — > — n + 2(z/ — 1). Similarly, the momentum of the state is given 
by 

/oo 2 
rfA (7(A) r 1 (A;/x) + E:g nQ (A°;/x) 
-CO i 



= P 9 r + Y,PnM), 
a=l 

where the breather momentum p n (X) is given by 

Pn(A) = - (-^n * (A) + g n (A ; /x) 
It is now easy to verify the important relation 

~PnW = e n (X) . 



(2.32) 



(2.33) 



(2.34) 



We remark that the following bootstrap-like relations are easily verified [1 1 



% % 
Ej(X+ -k) +s k (X - -j) = s j+h (X), 

s(X+ l -(v-l-j))+s(X- l -(v-l-j)) = ejiX). 



(2.35) 



Indeed, a hole (soliton) with rapidity A can be shown to have energy 7rsm/x s(A). We also 
remark that charge conjugation (C) and parity (V) eigenvalues can be readily computed using 
the methods described in Ref. ||14|| . Indeed, we find that the ground state is an eigenstate of 
C and V with eigenvalue (— ) 2". Moreover, an n-breather state has C = (— )"^ _n ; and if the 
rapidity is zero this state is also a parity eigenstate with V = (— )~2~ n . 

The preceding analysis, which is fairly standard, relied on only the first set ( 12.24] ) of 
Bethe Ansatz equations. In order to compute the two-breather S matrix, we also exploit the 



second set (|2.25 ) of Bethe Ansatz Eqs., which describes the centers of the breather strings. 
Forming the product iKLi an d taking the logarithm of both sides, we obtain 

&(A?) = J? , 
where h(X) is the new counting function 



(2.36) 



h(\) = -lN qni (\ ;t ,)-J2 

h=l 



£ r 2 (A^ - \ p ; /*) + £ Q2(X {niM ~ \ ft) 

13=1 h=l 



G2.37) 



and A^ 1 '^ = A + |(ni + 1 — 2Zi). We define the corresponding density <r(A) by 



We find 



(2.38) 



a ni (A ; ai) - £ / ^V(A') & 2 (A^) - A' ;//) + - £ a 2 (A^> - \t M ; //) 

h=i I/" 00 7V fa=i J 

1 711 f 2 "/3 712 

£ ni(A) + ^ £ E E (62 * ifi)(A^W - A^ } ) - £ a 2 (A^) - A^ ;/x) 

iV J X =1 [,0=1^=1 i 2 =l 



(2.39) 



In passing to the second line, we have used the result ( |2.27|) for cr(A). 



2.3 Breather bulk 5 matrix 

We define the two-breather S matrix S( ni ' n2 ' (X®, A 2 ) by the momentum quantization condi- 
tion 

(a?)* 5 (ni ' n2) (A?, A°) - l) |A°, A°) = , (2.40) 

where the breather momentum p n (X) is given by Eq. (|2.33j ). To compute the S matrix, we 
use the identity 



(2.41) 
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which immediately follows from Eqs. (|2.34| ) and (|2.38 ). Multiplying by i27iN, integrating 
with respect to A from — oo to A?, and noting the Bethe Ansatz Eq. (|2.36|) , we conclude that 
(up to a rapidity-independent phase factor) 

^ exp \i2nN J Xl dX (a- (A) - e m (A))| . (2.42) 

Substituting our result (|2.39|) for cx(A), we obtain 

ni ri2 

S (m,n 2 ) = p| JJ S (l,l) ^(ru,Ji) _ X (n 2 h^ 
h=lh=X 

ni n 2 , ■ v 

= II II [A -A + \{^-n 2 - 2h + 2l 2 )\ , (2.43) 

h=l h=l ' ~ 

where S^'^A) is given by 

cosh ({u- 3)f 



S^(X) 



cxp 



f duo , . cosh Uv - 3, 

1-2 — sinh (icuX) ) 

[ Jo u cosh Uv - V 



sinh cos (f(gj 

^h(^)-zcos(|(^f))' 



(2.44) 



This coincides with the sine- Gordon breather S matrix [II], 0, provided that we make the 
identification (3 2 = 8/i which we have already noted (|2.2|) . The breather S matrix has been 
obtained for the XXZ chain previously using the so-called physical Bethe Ansatz Eqs. in 

in. 



3 Open XXZ chain 

In this Section we consider the critical open XXZ spin chain with boundary magnetic fields 
h([i,£,±*) which are parallel to the symmetry axis Q 

w = 4 { £ («+i + «+i + cos ^ «+i) + h ^ > K + % , £+ K j , (3-1) 

where 

^(A* > = sin /i cot(//f) , (3.2) 

with < /i < | and e = ±1. For simplicity, we restrict to h(fi , £± ) < 0. 
3 The dependence on e is discussed in Appendix A. 
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Choosing again as the pseudovacuum the state with all spins up, the Bethe Ansatz 
equations are Q, []5f 

M 

e 2$ w_ 1 (A Q ; fi) e 2? ( E) _ 1 (A Q ; //) ei(A Q ; /i) 27V = J| e 2 (X a - Xp;fi) e 2 (X a + \p ; A*) , 

+ " 0=1 

a = l,---,M. (3.3) 

To streamline the notation, we shall often suppress the superscript (e) and thus write the 
boundary parameters as l;±. 

The energy is given by Eq. ( |2.3|) (plus terms that are independent of {A a }) and the S z 
eigenvalue is again given by Eq. ( f2.5| ). The requirement that Bethe Ansatz solutions corre- 
spond to independent Bethe Ansatz states leads to the restriction (see |§,|| and references 
therein) 

Re(X a )>0, A Q ^0,oo. (3.4) 

In addition to having the well-known "bulk" string solutions, the Bethe Ansatz Eqs. for 
the open chain also have "boundary" string solutions ||16|| . In particular, there are boundary 
1-strings A = ±z (t± — |) for | — | < £± < | (mod u). (See Appendix B.) For simplicity, 
we shall restrict £± so that such strings are absent, namely, 

- V -<i ±< l- V - (mod „) . (3.5) 

The lower bound comes from the restriction h([i,£±) < 0. 



3.1 One-breather state 

We consider again the attractive case e = —1. For values of £± in the range fl3.5| ), there are 
no boundary strings; and hence, the ground state is a sea of negative-parity 1-strings, as 
already discussed for the closed chain in Section 2.1. 

In order to compute the breather boundary S matrix, we consider the Bethe Ansatz state 
consisting of one breather Ag n ' Z ' ) in the sea. The corresponding Bethe Ansatz Eqs. read 

92^ + -i(X a ; aO 02£_-i(A Q ; n) ei(A Q ; fi) g\{K ; /i) 2Ar+1 

My n 
= - II e 2(A« - A/3 ; n) e 2 (X a + A^ ; //) J[ 92{K - Aq™' ; /i) g 2 (X a + Aq"' ; /i) , 

/3=1 Z=l 

a = 1 , • • • , Mf , (3.6) 
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i(Ao n '° ; y) e2c_-i(Ao n '° ; n) e 2 (2A[, n ' ; ft) e 1 (A^ ) ; /i) 



2iV 



A^i n 

II <MA^ - A, ; 2 (A< n ' l) + A, ; y) n e2 (Ai"'° - ; /*) e 2 (A^ + A?' fc) ; /*) 

0=1 



fe=l 



I = 1 , • • • , n . 



(3.7) 



The first set of equations ( |3.6| ) leads to the counting function 

h(X) = j(2iV + l)n(A ; fj) + qi (X ; n) + r 2€+ _x(A ; n) + r 2 £__ x (A ; M ) 

Aff n } 

£ fe(A - A/? ; /i) + g 2 (A + A, ; //)] - £ [r 2 (A - A^ ; /a) + r 2 (A + A< n '° ; »)] (3.8 



/3=1 



We define the corresponding density cr(A) as before ( 2.17 ). The restriction (3.4) on the Bethe 
Ansatz roots implies that we must pass from sums to integrals using §|,|| 



1 Ml ~ r°° 1 

^ £ <KA«) = I d\' a(X') g(X') - —g(0) 



(3.9) 



(plus terms that are of higher order in 1/N), where g(X) is an arbitrary function. We arrive 
in this way at the following linear integral equation for er(A): 

a(X) = -26i(A ;//)+/ dX' [a 2 (X - A' ; fjt) + a 2 (X + A' ; //)] tr(A') 

•/ o 



N 

1 n 
iV z=i 

Finally, defining the symmetric density cr s (A) by 

<x s (A) = 

we see that it is given by 



a x (A ; /i) + a 2 (A ; ^) + &i(A ; /i) + & 2€+ -i(A ; /i) + & 2€ __i(A ; fjt) 

A > 0. 



(3.10) 



a(A) A > 
ct(-A) A < 



(3.11) 



1 n 

a s {X) = 2 S (A)--E(^i(A-A^ ) ) + ir 1 (A + A^ ) ; 
iV i=i 

+ i- (R * (a x + a 2 + 6i + 6 2€+ _i + & 2 e— i)) (A) 



N 

where the Fourier transforms of R(X) and K n (X) are given by 

1 



R(lu) 



-1 + a 2 {u ; n) 
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(3.12) 



(3.13) 



and Eq. (|2.29|) , respectively. 

We turn now to the second set of Bethe Ansatz equations (|3.7|) . Forming the product 
Iir=i an d taking the logarithm of both sides, we obtain the counting function 



h(X) 



if n 

-\2Nq n {X ;/*) + £ [^ + -i(A (n '° ; /*) + g 2? _-i(A^ ; /*) + g 2 (2A^ ; /x) 

E E M A(n '° - A « > ^ + r *( A(n,i) + A ^ > a*) 



Z=l /3=1 



- E h(A (n '° - A (n ' fe) ; //) + g 2 (A w) + A (ri ' fc) ; //)] , (3.14) 



Lk=l 



where A*- n '^ = A + |(n + 1 — 21). We find that the corresponding density (7(A), defined as in 
Eq. ( |2.38| ), is given by 



i n ( 

a{X) = 2e n (X) + - E |& 2 (A (n '° ; //) + a^^A^ ; y) + a 2 ^(X^ ; y) + 2a 2 (2A^ ; /i) 

- (Kx * (en + a 2 + h + 6 2C+ _! + fr^-i)) (A^)| 

+1 E { (*i * 62) (A (n '° - 4^) + (#i * 62) (A (n '° + A^' fc) ) - 2a 2 (A^ + A^ fc ) ; y)\ 
^ i,k=i I ) 

(3.15) 

In obtaining this result, we have again used (|3.9| ) to pass from a sum to an integral, and 
then we have used our result ( p.!2j ) for the density u s (X) 

3.2 Breather boundary S matrix 

We define the boundary S matrix S^(X , £) for the breather Ajj n = A + |(n + 1 — 21) , 
I = 1 , . . . , n by the quantization condition 

(VWAo)iv sM( Ao) £_) SfW(Ao,e+) - l) |A > = 0, (3.16) 

where p n (X) is given by Eq. (|2.33|) . To compute the S matrix, we make use of the identity 



l^PnW + *(A) - 2e n (A) - = , (3.17) 
which is similar to Q2.41| ), and obtain (up to a rapidity- independent phase factor) 

S^ n \X , £-) S^(X , £+) ~ exp j^vriV jf*" dA (a(X) - 2e ft (A))| . (3.18) 
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Substituting our result (|3.15|) for <r(A), we obtain Q 

S {n) (\o,Z±) = Sf°(*>) S{ n \\o,Z±), (3.19) 



where 



71 IL 

4 n) (Ao) = II S£\\^ l) ) II S^(X^ l) + A< n ' fc) ) , (3.20) 



Kk 



and 



"Sq^C^o) — exp < — 2 f — sinh(2ic<jAo 

JO U) 



cosh (l) 


cosh (^f) 




cosh 


(("-1)1 


) cosh ((z/ 


-IV) 



sin (#T) -4(^ + 1) sin (^T) + 4(^1) - I ) sin ( 



27tAq 

2(^-1) 



(3 21) 

sin ( 2 1^ + 4(^ + f)sm(^- i ^ + f)sin(^ y + f) • 
We recall that S'( 1 ' 1 )(A) is the bulk two-breather S matrix ( |2.44 ), and that v = -. Moreover, 

n 

Sf ) (\o,t) = I[si 1 \\$ l > l \0, (3.22) 
i=i 



where 



00 dcu .^ u/r> . x cosh ((u - 2£)u) 
o) 



Si 1 ' (A , = exp I 2 / — sinh(2zo;A 
[jo c<j 

M»/-2m 
V 2(«/-l) J 



cosh ((z/ — l)u;) 



- sin J - cos 



sin(^f| - cos 



2(i/-l) 



(3.23) 



It can be shown that this result agrees with Ghoshal's bootstrap result |7|] for the breather 
boundary S matrix for the boundary sine-Gordon model with "fixed" boundary conditions, 
provided that we make the identification of bulk coupling constants (3 2 = 8// (see Eq. Q2.2p), 
as well as the identification of boundary parameters (] 

* = \ (2£ ( - 1} - p) ■ (3.24) 

In this formula we have restored the superscript (e) on the boundary parameter £. 

4 We assume that the terms which do not depend on £± contribute equally to S^ n ' (Ao , £-) an< l 

^ (,i) (Ao,e+). 

5 We denote the Ghoshal-Zamolodchikov (6),0 boundary parameter £ by x, in order to distinguish it from 
our boundary parameter £. We recall that Ghoshal-Zamolodchikov identify as "fixed" boundary conditions 
their case k — 0. Moreover, in the attractive case, their bulk coupling constant A = — 1 is related to our 
coupling constant v = — by A = z^ — 1; and their rapidity variable 9 is related to our variable Ao by 6 = Jzj;- 
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We remark that the appearance of the bulk S matrix in our expression for the boundary 
S matrix can be readily understood from the fact that an n-breather can be regarded as 
a bound state of n 1-breathers, which scatter among themselves upon reflection from the 
boundary This is illustrated in Figure p] for the case n — 2. A single line represents a 
1-breather, and so the 2-breather is represented by a double line. 




Figure 1: Scattering of a 2-breather from the boundary. 



3.3 Soliton boundary S matrix 

Although the main focus of this paper is on breather S matrices, we compute here the 
soliton boundary S matrix in order to further check the identification ( |3.24j ) of the boundary 
parameters. 

3.3.1 Attractive case (e = —1) 

We consider first the attractive case e = — 1, with £± in the range (|3.5| ), and so the ground 
state is a sea of negative-parity 1-strings. Following ||,|l(J, we consider the Bethe Ansatz 



state consisting of one hole with rapidity A in the sea, which has S z = +|. The counting 
function is 

h(X) = -^-|(2AT+l)r 1 (A;/i) + g 1 (A;/i)+r 2C+ _ 1 (A;/i)+r 25 __ 1 (A; / u) 

M- > 

- -\p;fj,)+ g 2 (A + L (3.25) 
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which leads to the density cr s (A) whose Fourier transform is given by [] 



C7» = 2s(uj)--J(co)(e wA + e- wA ) + 



N (-1 + a 2 (u ; fi)) 



+ d 2 (uj;n) + b 1 (uj;n) + b2 i+ -i(uj;n) + b2^_^i(uj;n) , (3.26) 
where J(u) is defined by 

l-a 2 (w;/i) 2sinh(f)cosh((z/-l)fJ 

We define the boundary S matrix S(X,£±) for the soliton by the quantization condition 

^ 2p( A)iV s{ ~ x j e _) ^ ; e+) _ i) | A) = , (3.28) 

where p(X) is given by 

p(A) = -(J*r 1 )(A)-r 1 (A;^). (3.29) 
The boundary S matrix has the diagonal form 



a(A,e±) 

/3(A,£±) 



(3.30) 



The matrix elements a(A , £) and /?(A , £) are the boundary scattering amplitudes for one-hole 
states with S z = +| and S z = — |, respectively. We compute these matrix elements with 
the help of the identity 



I> + «.(A)- 2s( A)-i>)=0. 



(3.31) 



We first compute a(A , £). We have (up to a rapidity-independent phase factor) 

a(A , £+) a(A , £_) ~ exp |i2vriV jf * (a s (A) - 2s(A)) dA j . (3.32) 
Substituting the result ( |3.26| ) for the root density and performing some algebra, we obtain 



a(A,0 ~ exp 



°° duJ ■ ufo- \\ rsmh(3(z/-l)f)sinh((z/-2)f 

— sinh I ziujX ) — r 

» w v y L sinh (f) sinh (2 (z/- l)w) 



sinh ((2f - l)w) 
2 sinh a; cosh ((z/ — l)u;) 



-cosh[7r(A + ^(2e-^))] 5 (A) 5i(A,0 

7T Z 



(3.33) 



6 Due to the presence of the hole, the prescription ( |3.9| ) for passing from sums to integrals has the additional 
term —jjg(X) on the right-hand-side. 
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where 



Sb(A) 



n 



r (-2iX + {v- l)(4n + 3) + 1) T (-2iA + - l)(4n + 1 



n=0 



T (2zA + (i/ - l)(4n + 3) + l) T (2zA + (i/ - l)(4n + 1) 
r (2iA + 4(i/ - l)n + l) T (2iA + 4(i/ - l)(n + 1 



T (-2iA + 4(i/ - l)n + l) T (-2iA + 4(i/ - l)(n + 1 



(3.34) 



and 

Si(A,£) 



n 



T (-iX + 2(i/ - l)n - ±(2£ - i/ - 1)) r (-iA + 2(i/ - l)n + §(2£ - v + 1] 



n=0 



X 



T (%\ + l)(2n + 2) - §(2£ - i/ - 1)) T (iA + (i/ - l)(2n + 2) + §(2£ - z/+ 1) 
r (iA + (i/ - l)(2n + 1) - ±(2£ -v-l))T (iX + (u - l)(2n + 1) + §(2£ - i/ + 1 



T (-iA + (i/ - l)(2n + 1) - §(2£ - v - 1)) r (-iA + (i/ - l)(2n + 1) + |(2^ - v + 1] 



(3.35) 



In order to compute /?(A,£), we must consider a one- hole state with S z = — ~. As 



explained in |]9fl,|[L0||, this can be achieved by working instead with the the pseudovacuum 
with all spins down, in which case the Bethe Ansatz Eqs. are given by (|3.3|) with the 
replacement £± — ► — £±. The corresponding density cr' s {ui) is given by Eq. ( |3.26| ) with the 
replacement 6 2 £ ± -i(^ ; /-0 — *• — &2£±+i(^ ; /•*)• We find 



o^(A) - cr s (A) = — (& 2 £_-*/(A ; tt) + 6 2 ? + - i ,(A ; tt)) 



(3.36) 



In obtaining this result, we have noted that for £± in the range (|3.5|) , the quantities &2£±-i (^ j A*) 
and &2£ ± +i(<^;A i ) are given by Eqs. ( |2.12|) and (|2.13|) , respectively. Since 



.7 c v 7? I v = ex P \ ^ N / (°s(A) - ff s (A)) dA 



(3.37) 



we conclude that 



cosh 



a(A,£) cosh 



7T 



(A - 1(2^ 



tt A + |(2£-i/ 



(3.38) 



The soliton boundary 5 matrix (|3.30|) , (|3.33|) - (|3.35|) ,( |3.38|) agrees with the bootstrap re- 
sult of Ghoshal-Zamolodchikov || for the boundary sine-Gordon model with "fixed" bound- 
ary conditions, provided the identification of boundary parameters ( 3.24 ) is again made. 
Fendley-Saleur M find a similar identification. 
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3.3.2 Repulsive case (e = +1) 



We consider finally the repulsive case e = +1, where there are no breathers. Here the ground 
state corresponds to a sea of positive parity 1-strings, i.e., real solutions of the Bethe Ansatz 
equations. For the Bethe Ansatz state with one hole of rapidity A in the sea, the counting 
function is 

h(X) = ^|(2Ar+l)g 1 (A;/i) + r 1 (A;/i) + g 25+ _ 1 (A;/i) + g 25 „_ 1 (A;/i) 



M~ s 

J2 - A^ ; n) + q 2 (X + Xp ; fi)\ \ 

8=1 ) 



(3.39) 



and the density a s (X) is given by 



+ a 2 (uj ; /j,)) 

+ a 2 (u ; fj) + bi(u ; (jl) + a 2 £ + -i(uj ; fi) + a 2€ __i(^ ; ji) 



(3.40) 



where now 



1 



S(UJ) 



J{u) 



aijuJj, /i) 

1 + a 2 (u ; fj,) 2 cosh (f) ' 

a 2 (u;n) sinh ((z/- 2) 



l + a 2 (^;/x) 2sinh((z/-l)f)cosh(f) 
Moreover, in the repulsive case, 

p(A) = -(J*fc)(A) + gi(A;^). 



(3.41) 



(3.42) 



Proceeding as in the attractive case, we find that the soliton boundary S matrix has the 
form (|3.30|) with matrix elements 



a 



(A , ~ exp <j 2 j — sinh [2iuj\j 



sinh(^) sinh ((i/- 2)1 
sinh (2a;) sinh — 1): 



+ 



sinh((z/-2£ + l)w) 



2 sinh ((i/ — l)u;) cosh a; 



I cosh[ ^_(A + -20)] 5 (A) S 1 (\,0 

7T V — 1 I 



(3.43) 



where 



fibw = n 



T (^(-2iA + 4n + 3) + l) r (^(-2iA + 4n + 1 



„=o I T (^(2zA + 4n + 3) + l) T (^(2iA + An + 1 
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r (+r(2i\ + An) + l)T (^r(2i\ + 4n + 4) . 
\ u ~ l - L \ u ~ l J— (3 44) 

T (^(-2?A + An) + l) T (^(-2iX + An + 4)) ' 
T (^(-iA + 2n - - 20) + |) T (^HA + 2n + \{y - 20) + \ 



S) I T (^(zA + 2n + 2 - \{v - 20) + \) V {^(iX + 2n + 2 + §(z/ - 20) + § 

r(^(»A + 2n + l-i(i/-2Q) + l)r(^(»A + 2n + l + l(i/-2e)) + |) 
T (iHA + 2n + 1 - \{v - 20) + \) T (A(-iA + 2n + 1 + \{y - 20) + \ 



(3.45) 



and 



cosh 


7T 




- 20); 


cosh 


7T 
V-\ 


(U 


K" - 20); 



(3.46) 



Comparing with Ghoshal-Zamolodchikov 0, we obtain the following identification of bound- 
ary parameters f\ 



x 



7T (u - 2^ +1 



2{y -1) 

The same identification was found in ||. We remark that, by setting 



(3.47) 



= 7T — /i . 



7/ 



(3.48) 



as in Appendix A, the formula ( |3.24p for x in the attractive case can be recast in the similar 
form 



x 



vr - 2Q) 
2(1/ - 1) 



(3.49) 



where z/ = 2L ! . 



4 Discussion 



We have formulated a systematic Bethe-Ansatz approach for computing breather S matrices 

for integrable quantum spin chains. We have used this approach to calculate the breather 

7 In the repulsive case, the Ghoshal-Zamolodchikov bulk coupling constant A is related to our coupling 
constant v by A = — ^y; and their rapidity variable 9 is related to our variable A by 6 = ir\. 
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boundary S matrix for the open XXZ spin chain with diagonal boundary fields. We have 
also directly computed the soliton boundary S matrix in the critical regime. 

Let us briefly compare our approach with that of other authors. Our approach is essen- 
tially a systematization of Korepin's analysis of the massive Thirring model. Key elements 
of our approach are the exploitation of the "second" set of Bethe Ansatz Eqs. (|2.25|) which 
describes the centers of the breather strings; and the use of the identity fl2.41| ). An analogous 
identity for holes was used by Andrei and Destri || to compute soliton S matrices. Fendley 
and Saleur || study boundary S matrices of the XXZ chain using an alternative approach 



based on the model's physical Bethe Ansatz equations [TT\. The identification of boundary 
5* matrices from the physical Bethe Ansatz Eqs. is not straightforward, especially in the 
repulsive case. Finally, we remark that the vertex-operator approach |1| has so far been 
restricted in applicability to the noncritical regime. 

While we have focused on the XXZ chain for simplicity, we expect that the same methods 
should be applicable to other models. Indeed, boundary S matrices for the critical 
open spin chain with diagonal boundary fields can be computed in this way [TO} . 
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A Dependence on e 

Following many authors (see, e.g., [j^|,flll||,|f[3||, ), we treat the full critical regime of the XXZ 
chain by restricting the anisotropy parameter /i to the range (0 , |), and introducing a new 
parameter e = ±1. We describe here this approach in detail, since there are some subtleties 
associated with it, such as the dependence on e in the expression ( |2.4| ) for the momentum 
and in the boundary parameters. 
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A.l Closed chain 



We take as our starting point the following definition of the critical XXZ closed chain Hamil- 
tonian 

1 N 

U = 7 E {«+! + + COS // - l) } , (A.l) 

4 n=l 

with < \j! < 7r. The repulsive regime corresponds to < // < |, while the attractive 
regime corresponds to f < // < 7r. The standard algebraic Bethe Ansatz procedure gives 

£ = -sinVE— T7 A , (A.2) 

^ ^ 1 cosh(2/x'A'J-cos/x'' V ; 

1« sinh/i'fA^ + f) 
P = + 7 E lQ g ■ r ; 777 (mod2vr), (A.3) 

* o=i smh fi' [\' a - 2) 

with 

/sinh/i / (A^ + f)\ iV Af sinh^(X' a -X'+i) 

/ { } =11 ^T 1 " (> a = l,---,M. (A.4) 

\sinh/x' ( A' Q -±W ^=1 sinh/x' ( A' a - X'p - i) 

These are simply the formulas of Section 2 for e = +1 with primes appended to /i and A Q . 
The principal observation is that the "duality" transformation 

jj, = 7T — // 

A. = $K + % (A.5) 



implies 



1 w 

« = tE {«+i + «+i " co S/U - l)} , (A.6) 

4 n=l 



M 



with 



£ = +sinVE — , , n ; x , ( A - 7 ) 

cosh(2/iA a ) — cos /i 

1 ^ sinh ^ f A a + |) 

p = 7rM- T ^log ) ^ (mod27r), (A.8) 

1 a =i sinh/i (A Q - | J 



sinh/x(A Q - §)/ smh/i(A Q - A^ - 
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The proof relies on elementary identities sinh //(A Q + | ) = — sinh /jf(\' a — |), sinh /i(A a — |) = 
sinh/i'(A^ + |), etc. 

The Bethe Ansatz Eqs. remain invariant under the transformation ( |A.5| ) for N even. 
Evidently, the attractive regime (-| < // < 7r) corresponds to < /i < |. The expressions 
( |A.7p ,( |A78D coincide with the corresponding formulas of Section 2 with e = — 1. Moreover, 
as is well-known []12|| , the Hamiltonian ( |A.6| ) with AT even can be mapped by a unitary 



transformation to the Hamiltonian (12.11) with e = — 1. 



A. 2 Open chain 

We define the critical XXZ open chain Hamiltonian by 

n = («+! + a n°n+i + cos /^' «+i) + sin^'cot(^)o-i + SU1 // COt (7 



n=l 



(A.10) 



with < // < 7r. The corresponding Bethe Ansatz Eqs. remain invariant under the transfor- 
mation (|A.5|) for any iV provided there is an accompanying transformation of the boundary 



parameters 



It follows that the Hamiltonian is equal to 

-JV-l 



1 f 1 

H = 4 \ Yl ( a n a n+l + a n a n+l ~ COS V «+l) ~ sin ^ Cot 04- ) ~ SHI /i COt(/i£ + ) (7^ \ , 

(A.12) 

which for any A" can be mapped by a unitary transformation to 

^ = -7 E ( a n a n+i + + cos p + sin /x cot(/x£_)of + sin fi cot ) a z N \ . 



, n=l 



(A.13) 

We conclude that the critical open chain can also be described by the Hamiltonian with 
< fi < | and e = ±1, where 

A+l) _ 

£± — ?± 

d" 1} = -(^-l)£i- (A.14) 
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B Boundary 1-strings 



The existence of boundary string solutions of the open-chain Bethe Ansatz equations was 
discussed in Ref. ||16fl . For completeness, we demonstrate here the existence of boundary 
1-strings following the approach used by Faddeev and Takhtajan |T7| to study bulk 2-strings. 
We therefore consider Eq. ( |3.3|) for the case M = 1 with iV — > oo, 



sinh a ( A + \ sinh a ( A + i(£ — ^ 

r %\ r — rr = 1 - 

sinh /if A — |W sinh/i f A — £(£ — |)J 

For simplicity, we have written boundary terms from just one boundary. Setting A = x + iy 
with x, y real, 

'sinh// (x + i(y + |))\ sinh/i (x + i(y + £ - |)) ^ ^ 
sinh /i fx + z(y — |)J y sinh /i + z(y — £ + |)) 

Multiplying by the complex conjugate and using the identity sinh(a + ib) sinh(a — ib) = 
sinh 2 a + sin 2 b, we obtain 

A 2N B = 1, (B.3) 

where 

sinh 2 fix + sin 2 fi(y + §) sinh 2 /ix + sin 2 /i(y + £ — ~) 

= o 2 — 7 TT ' = 2 2 — ~i TT • (-^ ,z v 

sinh fix + sin fi(y — ^) sinh fix + sin — £ + 

Evidently, there is a periodicity y y + E . We therefore consider two cases: 
• Case I: 

< y < — (mod -) (B.5) 
2fi fi 

It follows that sin 2 fi{y + ~) > sin 2 /i(y — |); therefore A > 1, and hence A 2N — > oo for 
iV — > oo. Eq. QB.3] ) then implies B — > 0. That is, for iV — > oo, 

x = 0, y = -(£_!) (mod-). (B.6) 

The restriction ( |B.5| ) then implies 

-- — <£<- mod-. B.7 
2 2/i 2 /i 
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Case II: 

- — < y < (mod -) (B.8) 
2/i n 

Then sin 2 ^1(2/ + ~) < sin 2 yu(?/ — ~); therefore A < 1, and hence A 2Ar — > for N — > oo. 
Eq. (|B.3|) then implies 5 — ► oo. That is, for AT — > oo, 

x = 0, y = (mod-). (B.9) 



The restriction ( p.8| ) leads again to the condition (|B . 7| j . 



In conclusion, the Bethe Ansatz Eqs. have the boundary 1-string solutions A = ±i (£ — | 
when £ satisfies the condition (p3.7| ). Boundary strings of longer length are also studied in 
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